function [cube] = radiance2reflectance(cube,bands,whitecube,whitebands,whitemask,white_reflectance_filename)
%get_reflectance 此处显示有关此函数的摘要
%   radiace
%   white_radiance_filename

% wr = csvread(white_radiance_filename,53);
% wr(:,3) = [];
[~,index,~] = intersect(whitebands,bands);
whitecube = whitecube(:,:,index);

wf = csvread(white_reflectance_filename,1);
wf(:,3) = [];
[~,index,~] = intersect(wf(:,1),bands);
wf = wf(index,2);
for i=1:size(cube,3)
    cube(:,:,i) = whitemask.*cube(:,:,i)./(whitecube(:,:,i)+eps)/wf(i)*100;
end

end

